Efficient Decoupling Schemes Based on Hamilton Cycles 
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Decoupling tlie interactions in a spin network governed by a pair-interaction Hamiltonian is a 
well-studied problem. Combinatorial schemes for decoupling and for manipulating the couplings 
of Hamiltonians have been developed which use selective pulses. In this paper we consider an 
additional requirement on these pulse sequences: as few different control operations as possible 
should be used. This requirement is motivated by the fact that optimizing each individual selective 
pulse will be expensive, i. e., it is desirable to use as few different selective pulses as possible. For an 
arbitrary d-dimensional system we show that the ability to implement only two control operations 
is sufficient to turn off the time evolution. In case of a bipartite system with local control we show 
that four different control operations are sufficient. Turning to networks consisting of several d- 
dimensional nodes which are governed by a pair-interaction Hamiltonian, we show that decoupling 
can be achieved if one is able to control a number of different control operations which is logarithmic 
in the number of nodes. 

PACS numbers: 03.67.Lx, 03.65.Fd, 03.67.-a 
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I. INTRODUCTION 

The problem of simulating pair-interaction Hamilto- 
nians has been studied intensively Q, S S S S 01 ■ It 
has been shown that starting from a given entangling 
bipartite Hamiltonian any bipartite Hamiltonian can be 
simulated . This can be extended to networks of qubits 
0, 0, and even to networks of higher dimensional sys- 
tems ^ m . The underlying idea is to control the system 
by applying hard pulses to the individual nodes of the 
network. This parallels the methods developed for uni- 
versal control of open quantum systems 0, B 0, ■ De- 
coupling of arbitrary interactions is an important prim- 
itive for the simulation. This can be achieved by using 
combinatorial constructions such as triples of Hadamard 
matrices '3| and orthogonal arrays The control op- 
erations in these schemes are local and have to be applied 
selectively to the nodes. In the following we give a brief 
account of the general setting. 

We assume that i7 is a pair-interaction Hamiltonian 
of a system which consists of N coupled spin-i particles. 

Let aa\ where a € {x, y, z} and i € {1, . . . , N} denote 
the Pauli operators acting locally on qubit i. We can 
write H in the form 
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with indices a, P € {x, y, z}. Hence the coupling strength 
between the different qubits in this iV-spin network is 
given by J^'^. The problem at hand is to decouple inter- 
actions in a Hamiltonian of the form cq. or more gen- 
erally to simulate another Hamiltonian H by the given 



one. The framework for these simulations 
Hamiltonian theory [T^ . 

A special case is given by a system which has azaz 
interactions only. Here selective ax pulses applied to the 
nodes are sufficient to turn off the interactions and hence 
to simulate any desired Hamiltonian P, Q . 

For general pair-interaction Hamiltonians the decou- 
pling problem is harder. Nevertheless, in it was 
shown how to achieve decoupling using orthogonal ar- 
rays which are objects studied in combinatorial theory. 
The approach taken in |lll[generalizes schemes obtained 
from Hamadard matrices |2 ■ AH these methods can 
be thought of as generalizations of well-known techniques 
for decoupling and refocusing used in nuclear-magnetic- 
resonance theory [13 . IT3 | . In higher-dimensional systems 
decoupling can be achieved by means of selective pulses 
which are derived from orthogonal arrays 0, ^Jl- The 
selective pulses employed in these schemes will be quite 
demanding in experimental realizations. Hence it would 
be desirable to have as few different selective pulses as 
possible in order to minimize the optimization overhead 
necessary to implement each one of them. This motivates 
the question to search for decoupling schemes in which 
the pulses can be arranged in such a way that only few 
different pulses are used. 

After a brief introduction into the framework we will 
consider schemes for one d-dimensional node which will 
work with two different pulses, schemes for a bipartite 
system of two c?-dimensional nodes which will work with 
four different pulses, and finally schemes for networks of 
n nodes which are d-dimensional. The latter schemes 
make use of 0(log7i) different pulses. 



II. DECOUPLING SCHEMES 
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The underlying model of the decoupling schemes de- 
scribed in this paper is average Hamiltonian theory [l3j |. 
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the necessary parts of which we briefly describe next. As- 
sume that the Hamihonian H is of the form . We can 
apply the sequence 

e-^r^HVn...e-'^^"V2e-''-^"Vi, (2) 

with relative times G R and local unitaries Vi for 
i = l,...,n. Note that the requirement rir=i ~ -'- 
is necessary in order for average Hamiltonian theory to 
hold and to determine the terms in the Magnus expansion 
|14j of the piece- wise constant evolution given in eq. l(2Jl • 
We can rewrite ^ in the form 

(K . . . Vi)^e~"-" . . . {V2Vi)^e-"^"{V2Vi) Vle-'^'"Vu 

i. e., as a product in which all factors are given by Hamil- 
tonians Hj — UjHUj, where Uj := 01=1 ^fe (ordered 
from the right). Hence, the time evolution is divided into 
n intervals in each of which we have a conjugated time 
evolution Uje~'^'^''^Uj with respect to the basis given by 
Uj, i. e., this is the "toggling- frame" form of the sequence 
p^. We can still rewrite this using Uj = 111=1 and 
obtain the new sequence 

Ule-'^""UnUl^, . . . ule-'^'"U2Ule~'^'"Ui. (3) 

This is the form of the pulse sequence we actually work 
with in the following. The effective Hamiltonian corre- 
sponding to eq. i.e., the first term in the Magnus 
expansion, is given by ^"^1 ujHUi. If the unitary oper- 
ators Ui applied in the toggling-frame have the property 
that for each pair (fc, I) of nodes the respective local oper- 
ations applied to spins k and I run through the elements 
of a unitary operator basis B 0,0,0] of the subsystem 
given by k and I, then we obtain that the time evolution 
of the system is stopped. Indeed, on the subsystem we 
then obtain the average Hamiltonian ^ 

which is zero since ti{H) = 0. Unitary operator bases 
B exist in any dimension. Indeed, let d > 2 be the di- 
mension of the system, we explicitly define a basis for 
the vector space of C'^^'' matrices: let cud = exp(27ri/c?), 

<^x = J2i=o N) (« + = EjJo N) where aU in- 

dices are computed modulo d. Then B = {(r^cri : i,j — 
0; . . . , d — 1} is an operator basis which will be referred 
to as the Pauli basis. 

We continue with an observation on the sequence 
while the matrices Ui have to run through the list of all 
elements of the operator basis B, the list of quotients 
Ui+iU} can be considerably smaller. The following ex- 
ample shows that for a system consisting of two d di- 
mensional subsystems it is enough to be able to apply 
four different operators in order to switch off the time 
evolution. 

Example. In Figure ^ we have shown how to label 
the elements of the Pauli basis for a d dimensional sys- 
tem in such a way that we reach all elements by just 
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FIG. 1: Decoupling of a Hamiltonian on can be achieved 
using two different pulses (the case d = 6 is shown). Pulses 
correspond to elements in Z| which is generated up to phases 
by Gx and G:,. The left square shows the sequence in which 
the pulses have to be applied. In each square a horizontal 
move corresponds to multiplication with o-i'' and a vertical 
move I to multiplication with cri'' for i = 1,2. An alternative 
sequence is given in the right square. Decoupling of a bipartite 
Hamiltonian on C* (g) can be achieved using four different 
pulses. Here pulses correspond to elements in = x Z^ 
generated by cy'i\ cri*' for i = 1,2. A cyclic sequence can 
be obtained by performing one step in the first copy of Z^ 
followed by one full cycle in the second copy and so on. 



multiplying one generator at a time. It is also shown 
how the elements of an operator basis for a bipartite 
system can be arranged. This corresponds to an enu- 
meration of the elements of = {S)/{{S) n C*), where 

S := {alc^^ , cri^'' , tri^^ , ai^''}. This enumeration can be cast 
in terms of Hamilton cycles in Cayley graphs (for the ba- 
sic notions of graph theory we refer to [13 )■ Let G be a 
group and let 5" be a set of generators of G. Then the 
Cayley graph F(G', S) — {V, E) is the directed graph with 
vertices V = {vg : g e G} labeled by the group elements. 
The edges E are as follows: there is a directed edge from 
Vg to Vh if there exists s G S such that h = sg. Each 
Cayley graph is regular of degree \S\ and its structure 
depends on the specific choice of the generating set S. 

A tour through the vertices of a graph which visits 
each vertex precisely once is called a Hamilton path. In 
case there is an edge between the end node and the start 
node we can close the cycle and obtain a Hamilton cycle. 
The sequences shown in Figure ^ define Hamilton cycles 
in the Cayley graph T{Z%S). 

While for graphs like the Cayley graph of Z^^ Hamilton 
cycles can easily constructed directly, the general ques- 
tion whether Hamilton cycles exist in arbitrary Cayley 
graphs is open 0|. In general the problem to decide 
whether a given graph contains a Hamilton cycle is a 
difficult problem and of interest in computer science and 
optimization. It is one of the classical NP complete prob- 
lems ,2Qj, i.e., it is believed that no polynomial time ex- 
ists for this problem. 

For the rest of the paper we use the fact that there 
are Hamilton cycles in the Cayley graphs r(G, S) for the 
abelian groups G — Z^" and the set S of generators is 
given by the 2n coordinate vectors. These can be ob- 
tained by generalizing the construction given in Fig. ^ 
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III. ORTHOGONAL ARRAYS 

In the design of statistical experiments which de- 
pend on several factors one is often forced to work with 
an incomplete subset of the set of all possible combina- 
tions of factors. Orthogonal arrays provide a way to plan 
such experiments systematically. An orthogonal array is 
a (in general non-square) matrix with entries from a finite 
set S of symbols. The properties of an array are deter- 
mined by a set of characteristic parameters for which the 
notation OA{N, n, s, t) is used |23|. An n x N matrix A 
with entries from the set S is an orthogonal array with 
s := \S\ levels and strength t iff every tx N subarray of A 
contains each possible i-tuple of elements in 1 5*1 precisely 
once as a column of A. The shorthand notation for this 
is that A is an OA{N, n, s, t). 

The rows of this matrix are also called factors and 
correspond to nodes of the network. The columns are 
also called runs and correspond to the time slots of the 
scheme, i. e., to the pulses to be applied. The most impor- 
tant parameter is the strength t of the array. In the con- 
text of simulation of pair-interaction Hamiltonians the 
strength has to be i = 2. 

In it has been shown that given an orthogonal 
array OA[N, n, 4, 2), an arbitrary and possibly unknown 
Hamiltonian H, which describes the pair-interactions in 
an n qubit network, can be decoupled using TV pulses. 

We now show that in this case the pulse sequence can 
be arranged in such a way that only a minimal number 
of different pulses have to be applied. First, we briefly 
recall some basic facts about error correcting codes |23| 
since they will feature in the subsequent constructions of 
orthogonal arrays. A linear code over the flnite field 
is a /c-dimensional subspace of the vector space F^. The 
metric on the space F^ is called the Hamming weight 
which for x = (xi, . . . ,a;„) G F^ is defined by wt(a;) := 
\{i G {1, . . . , n} : Xi ^ 0}|. The minimum distance of a 
linear code C is defined by c? = dmin ■— min{wt(c) : c G 
C}. As a shorthand we abbreviate this by saying that C 
is an [n, k, d]q code. We need one more definition which 
is the dual code C-^ of C defined by C-^ := {a; G F^ : 
a; • 2/ = for aU ?/ G F^}. 

The following theorem j2^ Theorem 4.6] establishes 
a connection between orthogonal arrays and error cor- 
recting codes. In fact this is one of the most prolific 
construction for orthogonal arrays known. 

Theorem 1 Let C be a linear [n,k,d]q code over Fq. 
Let d^ be the minimum distance of C^. Arrange the 
codewords of C into the columns of a matrix A G F^'^'^ . 
Then A is an OA{q'' , n, q, d-^ — 1). 

Each column of the orthogonal array corresponds to 
one of the pulses Ui in eq. (O and we want to keep the 
set {Ui+iUj} as small as possible. We are going to show 
next that for the arrays constructed from Theorem^this 
is always possible. 



IV. DECOUPLING USING GRAY CODES 

We illustrate Theorem ^ by a scheme for a small net- 
work consisting of 5 qubits. The flnite field needed for 
this case is the field F4 — {0,l,uj,uJ} of four elements in 
which the relation 1 -|- + cJ = holds. The following is 
a generator matrix for the quadratic residue code (23| C 
over F4 with parameters [5,2,4]. 

" 1 1 w w" 
1 IZJ tZJ 1 

The dual code of this code again is defined over 
F4 and has parameters [5,3,3]. By taking all 4^ — 16 
code-words in C as the columns of a matrix we obtain 
an OA(16, 5, 4, 2). We now address the question how to 
arrange the columns of this 4 x 16 matrix. First, recall 
the Gray codes [23.l25l | are Hamilton cycles for the groups 
F2 where the generator sets are the coordinate functions. 
For F2 a particular Gray code is given by the following 
sequence of binary strings: 

■0 00000001111111 1" 
0000111111110000 , . 

0011110000111100 ■ ^ ' 

. 1 1 1 1 1 1 1 1 . 

Observe that from column q to q+i exactly one position 
is fiippcd. Overall we obtain a cyclic sequence which 
runs through all 16 elements of Fj exactly once. The 
order defined by the Gray code also defines an order of 
the vectors in F4. To obtain the corresponding elements 
in F4 we make the identification (0,0) 1— > 0, (0,1) 1— > 1, 
(1,0) I— > uj, and (1, 1) 1— > tZJ. Using this identification we 
can now map all elements of F4 to elements of F4 using 
the linear code C. We obtain the following scheme for a 
system of five qubits: 

"00 1 1 1 IujUUujujujlouj' 

1 u lo lo lu 1 i u lo lo lu 1 

u lu 1 lu lu 1 uj 1 lu lu 1 lu 

u lu 1 lu 1 lu lu 1 lu u lu 1 

_ 1 u lu 1 lu lu uj u 1 lu lu 1. 

Hence there are four basic pulses tti , . . . , 7r4 correspond- 
ing to those columns in Q which are elementary vec- 
tors ei,...,e4, i.e., columns 2, 4, 8, and 16. The cor- 
responding operators are tti = I2 ® ci (X> us ^ us (8) ti, 

TT2 = I2 "Xl 0-2 <X) CTi (X) tJl <X) 0-2 , 7^3 = CTi (g) I2 I?) CTi (g) (73 (g) (T3 , 

7r4 = (72 <8) I2 €5 172 ® (Ti (8) <Ti . The order in which they 
have to be applied can be read off from the transitions 
in Q), i. e., we obtain the following pulse sequence: 

TI"!, 71'2, TTl, TTs, TTl , 772 , TTl , 7r4 , TI"! , 7r2 , TTl , TTa , TTl , 772 , TTl , 774 . 

We now turn to the question how to generalize this 
to more general networks and show that pulse sequences 
obtained from linear codes have the property that the 
number of different pulses grows logarithmically with the 
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number of spins. In order to apply Theorem^we have to 
find linear codes over F4 for which the minimum distance 
of the dual code is at least 3. 

Let g be a prime power and let m G N. The 
Hamming code Hq „i of length n := (q™ — 1) is 
a single-error correcting linear code with parameters 
[n,n — m, 3]g. The corresponding dual code -ff^^ has 
parameters [n,m,q"^~^]. By specializing q = 4 and 
by using Theorem ^ for -ff|^„j we therefore obtain or- 
thogonal arrays with parameters OA(iV, n, 4, 2), where 
n = (4™-l)/3 and N = 4"-™ for any choice of m G N. 

The procedure to obtain a pulse sequence for a network 
of n spin-i particles, where no is an arbitrary natural 
number, i. e., not necessarily of the form n = (4™ — l)/3 
is as follows: first let m G N be the unique integer 
such that no < ^—3—^ < 4no. Then construct the or- 
thogonal array with parameters 0^(4™, (4™ — 1)/3, 4, 2). 
The columns of this orthogonal array are codewords of 
^4^m Q F4'* We will use only the first no rows of 

this ^-3—^ — m X ^-3^ array. Since H^^^ = FJ*, we can 
find the desired Hamilton cycle by choosing a Gray code 
on FI™, i.e., we can decouple using 2m different pulses. 

Higher Dimensional Systems. — So far we have only 
considered spin-i particles, i.e., networks consisting of 
qubits. Methods exist to decouple Hamiltonians also in 
case the dimension of the individual nodes is greater than 
two |3) ^ . Basically, for one node of dimension d> 2 the 
requirement to switch off the time evolution of this node 
is to apply all d? elements of a unitary operator basis 
{Uij : i,j = 0,...,d—l}. For any pair of nodes we 
have to apply all d"* elements of a tensor product basis 
{Uij Uk,e ■ i,j,k,£ = 0,...,d—l}. These pulses can 



be arranged like in the previous schemes using orthogonal 
arrays. By combining the construction of orthogonal ar- 
rays from Hamming codes over the alphabet F22Q , where 
a G N, together with the Gray codes for F22m, which 
can be constructed as in Figure ^ we obtain schemes for 
networks of n = (2^™ — l)/3 nodes (each of dimension 
2") which use only 41ogn different pulses. 

Conclusions and Discussion. — We have shown that 
pulse sequences for decoupling in networks of qubits and 
higher-dimensional systems can be arranged in such a 
way that the number of different pulses needed is small. 
Using orthogonal arrays derived from Hamming codes 
over finite fields it has been shown that the number of 
different operations can be chosen to be of the order 
0(log n), where n is the number of nodes in the network. 
The method presented in this paper is based on Hamilton 
cycles in the Cayley graph of the group defined by the 
columns of an orthogonal array. The related graph theo- 
retical concept of Eulerian cycles has been used recently 
for the problem of designing robust schemes for decoher- 
ence control [26L |2^ . Since in this case the Hamiltonian 
of the system can be arbitrary, i.e., not necessarily of 
pair-interaction type, the task to decouple from an envi- 
ronment leads to a different control scenario. It would be 
of interest to determine whether the methods described 
in this paper can also be used to derive efficient schemes 
for robust Eulerian dynamical coupling. 
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